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QUASI-SLENDER BODY THEORY FOR SLOWLY OSCILLATING 
BODIES OF REVOLUTION IN SUPERSONIC FLOW 


SUM/VIARY 


An analysis is presented which accounts for body shape and Mach number 
dependence of the aerodynamic forces on slowly oscillating pointed bodies of revo- 
lution in supersonic flow. By using body -fixed coordinates, the first order velocity 
potential is expanded in an elementary fashion for small radial distances from the 
body. This expansion is equivalent to the Adams-Sears expansion technique using 
Laplace or Fourier transforms. Analytical closed-form expressions for the sta- 
bility derivatives are obtained. These results are valid to second-order in terms 
of Mach number and body thickness ratio. 

Three specific body shapes are analyzed in detail: a cone, a convex 
parabolic ogive, and a concave parabolic ogive. Numerical results for the four 
stability derivatives , C + , and C + C are presented to 

Q! q; Q! q o; q 

illustrate their dependence on Mach number and body thickness ratio. Compari- 
sons with other theories are given, and the range of validity of the present results 
is discussed. Comments on application of quasi -slender -body theory to other body 
shapes are made. 


INTRODUCTION 


For the calculation of static and dynamic stability characteristics of large 
rocket vehicles in the high dynamic pressure portion of the ascent, which is usually 
the low supersonic flow regime, there exists a need for an easy-to-apply theoreti- 
cal method that will show explicitly the effects of Mach number and body geometry. 
Munk*s well-known slender -body concept [1] and its generalization to unsteady flow 
[2] are often found to be too inaccurate. Consequently, in recent years many in- 
vestigators have tried to develop more accurate methods by either seeking a solu- 
tion to the full linearized potential equation [3-11] or by attempting to solve the 
nonlinear equation [12-16], 

Extensions of the Karman -Moore method [17] to the case of slowly os- 
cillating bodies of revolution are given in References 8, 9 and 10. A similar ex- 
tension of the Oswatitsch-Erdmann procedure [18] is proposed in Reference 11. 
References 12, 13 and 16 consider linearized time -dependent perturbations about 
the nonlinear axisymmetric flow. The resulting equations are then solved by the 


method of characteristics. All of these approaches lead to purely numerical pro- 
cedures and require a substantial amount of preparatory work. 

A quite different nonlinear approach was developed by Revell f 14] which is 
an extension of LighthUl's second -order slender -body theory fl9] to slowly oscil- 
lating bodies of revolution. This method makes possible the determination of 
analytical closed -form expressions for stability derivatives. However, the eval- 
uation becomes rather unwieldy for bodies with arbitrary meridian profiles. 

This report presents a linearized approach which allows the derivation of 
closed-form expressions for the stability derivatives. By using body -fixed coor- 
dinates, a solution to the first-order potential equation is sought. For this pur- 
pose the velocity potential, as given by Dor ranee [3] , is expanded in an elemen- 
tary fashion for small radial distances from the body. Such an expansion is 
equivalent to the iterative method of Adams -Sears f 20] so that the doublet strength 
may be obtained in closed form. This procedure provides a consistent generali- 
zation over slender 4body theory, which is the first term in the expansion and 
hence prompts the name ’’quasi -slender body theory. ” The higher order terms 
represent the dependence of the solution on Mach number and body geometry. 


The present work considers only the next higher order term beyond 
slender-body theory, i. e. , only two terms in the Adams -Sears iteration. One 
of the chief aims in the derivation given here is to achieve ’’consistency" in the 
sense that all so-called second -order terms which arise from the velocity po- 
tential and pressure coefficient will be included in the stability derivatives. 

To date, for bodies of revolution, the Adams-Sears method has been 
applied to only one configuration, the slowly oscillating convex parabolic ogive 
body, in transonic flow by Landahl [6] and supersonic flow by Zartarian-Ashley 
[7] . In both of these treatments the aerodynamic forces were obtained by 
’linearized" momentum considerations. 

The expressions for stability derivatives which are worked out in this 
report are applicable to smooth slender pointed -nose bodies of revolution with 
arbitrary meridian profile. Three numerical examples are given to illustrate 
the effects of Mach number, differing body geometry and thickness ratio. 

These examples are a cone, a convex parabolic ogive, and a concave parabolic 
ogive. Several comparisons of the present results with other theories are pre- 
sented also. 

In the following treatment, dimensionless variables are used with dis- 
tances referred to body length, f * , velocities to free stream speed, U* , 
velocity potentials to product of body length with free stream speed, and time 
to body length divided by free stream speed. Starred quantities are dimen- 
sional, whereas unstarred quantities are dimensionless. 


PROBLEM FORMULATION 


Potential Equation 


Consider a harmonically pitching pointed body of revolution which is set 
in a steady uniform supersonic stream. Both the amplitude and frequency of 
oscillation are assumed to be small. A body -fixed cylindrical coordinate sys- 
tem, as shown in Figure 1, is used to describe the problem. The body is re- 
quired to be sufficiently slender so that shock wave effects are negligible. 
Therefore, a velocity potential (x, r, 0, t) exists. Following Revell [14] 
a body -fixed perturbation potential ^ (x, r, 0, t) may be defined by 


(x, r , 0,t) = (x-a) cos 6 + r sin 6 cos 0 (x,r, 0,t) , 

where 

5 = 6fl cos k t = 

The velocity components in the (x, r, 6) directions are then given by 


( 1 ) 


u = = cos 6 + ^ 

X X 

V = = sin 6 cos 0 ^ 

r r 

1 ^ 1 
w = — S2 = - sin 6 sin 0 + — ^ ^ 

r 0 r 0 


( 2 ) 


Nonlinear terms in the potential equation will be neglected in the following 
treatment. With this restriction ^ is considered to be governed by the first-order 
unsteady equation. 


^ - $ 
XX rr 


I 

V 


r' 


- - * - -2 ^ee <2*^,+ V =“• 


(3) 


For steady flow the superiority of the first-order equation over the linearized 
equation has been established by Van Dyke [21], For unsteady flow, to the 
authors’ knowledge, this point has not been established so that the superiority 
of equation (3) must be taken as tentative. 

In this treatment body -fixed axes have been chosen because the tangency 
condition and pressure coefficient are less cumbersome to apply than in a wind- 
fixed coordinate system. The wind-fixed ^stem is complicated by the zero angle- 
of -attack solution entering into the cross-flow boundary conditions through the 
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Taylor series transferral of the tangency condition from the pitching body to the 
mean body [22] . In a body -fixed system this complication is absent. 

For small amplitude and low frequencies, following C. Heinz [23] , the 
perturbation potential may be written 

$(x,r, 0,t) = 0(x,r) + [ 6 • (x,r) + 6 * A (x,r) ] cos d, (4) 

where 6 = (ik6) jj, <{> is the steady -state zero angle -of -attack perturbation po- 
tential, and ^ and X are the in-phase and out -of -phase cross-flow perturbation 
potentials, respectively. 


Boundary Conditions 

Upstream of the body all perturbations vanish; hence, 
(p = !p = X = 0 for X - /3r ^ 0 

p = $ = X =0 for X - /3 r = 0 

XXX 


(5) 

( 6 ) 


In addition, the flow must be tangent to the body, which in body -fixed coordi- 
nates may be expressed as 


p = R' (x) • + p } 

r X 


> at r = R(x) . 


\ = - [ (x-a) + R(x) R’(x) ] +R’(x) • j 


(7) 


These are the so-called exact tangency conditions for a non-deforming body. 


Normal Force and Pitching Moment 

The normal force and pitching moment, referred to the body base area 
Q=!' (£'•' ) = TT i''' ^ and body length , are given by the relations 

in 

C = - ''t: f f RC (x,R,0,t) cos 0 d0 dx (8) 

N Q(l) n P 
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and 


C_ = 


1 7T 


M Q (1) 


J J (x-a) RC^(x,R,0,t) cos 0 d0 


dx 


0 0 


1 7T 


Q (1) 


f J* R^ R’ C (x,R, 0,t) cos 0 d0 dx, 


(9) 


0 0 


where the pitching moment is taken about an axis at x = a. The last integral in 
the moment coefficient is the contribution of forces parallel to the body axis. 

Analogous to the expansion of the perturbation potential the pressure 
coefficient may be expanded for small frequency and amplitude : 


Cp(x,R,0,t) = Cp (x,R) + [6 Cp (x.R) + 6 (x,R)] cos 0 . 


(10) 


Then the normal force and pitching moment reduce to 


7T 


~ ^ Q(l) 


C =6 
M 


+ 6 


Q (1) 


Q 


f R C (x,R) dx - 6 


P ' • Q(l) 


f R C (x,R) dx (11) 

n 


0 P 2 


X 

f (x-a) RC^ (x,R)dx+ / R^ R ' (x,R) dx 


0 

1 


f (x-a)RC (x,R) dx + prfr- / R^ R’C (x,R)dx .(12) 

(1) -"o Pg y(l) P 2 J 


The flutter derivatives and are just the coefficients 

^ JN^ IN* • 

of 6 and 6 tn equations ( 11) and ( 12) . 


Expansion of the Pressure Coefficient 


The exact isentropic pressure coefficient relation is derived from 
Bernoulli's equation for unsteady flow [14, 24] : 


p y M 


^ G 


- 1 


(13) 
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where 


G 


1 + r 1 - ( + w^ + 2 ) ] 

2 I ' j' j 


(14) 


and subscript t denotes partial differentiation with respect to time in a wind- 
fixed coordinate system. Revell [14] has shown that for small angles of pitch 
the partial derivative with respect to t may be expressed in terms of body- 
fixed variables by 


8 
8 T 


8 
8 t 


- 6 r cos 6 


8x 


+ 6 


(x-a) 


r 8 

sin 6 

8 

[cos e - 

r 

80 _ 


. (15) 


By using equation (15) together with equation (2) , which relates the velocity com- 
ponents to the perturbation potential, and equation (4) for the expansion of the 
perturbation potential for small amplitude and low frequencies, G becomes 

G = 1 - ^ [Fo + 2 ( 6 Fj + 6 Fa ) cos e ] , ( 16) 

where 


Fa = 2 +0^ +0^ 

^ X X r 


Fj = (1 + <t>^) + (1 + 


Fa = X {i + <j> )+il)+(l> {X + x-a) -V<P 
XX r r X 


The pressure coefficient may be expanded using the binomial theorem. To 
first-order in 6 and 6, this gives 


C 

P 


-Fo (1 - — Fo) - 2 (1 - — Fo )( 6 • Fi + 6 . Fa) cos 6 . 


(17) 


We now put the exact tangency conditions, equation (7) , into the definitions of 
Fi and Fa, then substitute into equation (17) and retain terms up through second- 
order in Mach number and body thickness. By comparison of this result with 
equation (10) , we obtain the following expressions for the in-phase and out-of- 
phase lifting pressure coefficients evaluated at the body surface r = R(x) . 


C (x,R) = -2 il) 
p ^x 

1 


i - + a - R 


.2 


(18) 


6 



( 19 ) 


C (x,R) = 


2 

+ ,p a <P^~ ^ R’2) - RR»2 -Rct>^ 


These expressions are consistent with the order to which the potential equation 
will be solved. 

SOLUTION OF THE POTENTIAL EQUATION 


Assuming harmonic time dependence, 

$ (x,r,0,t) = <#>(x,r) + <^>(x,r,0) • e^^* 
the cross -flow potential $(x,r,6) satisfies 

-A A IAIA a a, a 

fl2 (P - $ _ 1^ 0 - -2 + m 2 (2 i k 0 - r <#> ) = 0 

^ ^xx ^rr r r r‘‘ ee 


X 


( 20 ) 


( 21 ) 


The well known dipole solution [25] is 

cos 6 of -iju(x~4) cos K p 


4> (x,r,0) = - 
dipole 

where 

P 


27t Or 


’ e-"' 


= 's/(x-|)2- ;S2 j,2 


K = 


M = 


k M 
k m2 


(22) 


B5t distributing such oscillating dipoles along the x-axis, we obtain the 
general solution for an oscillating body of revolution. 


(j) (x,r,0) = - 


cos 6 8 

27t 0 r 


x-/3r 


J ^ F(f ) cos K p jJL {x- ^ ^ ^ 


(23) 
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where F( is the dipole distribution which must be determined from the boundary 
conditions at the body surface. 


For a slowly oscillating body of revolution, we may expand this solution 
with respect to the frequency and retain only terms up through the first power of 
the frequency; thus, 


A 

<p (x,r,6) 


cos 9 
2 7T 


9 
8 r 



0 


F(g) ri -iM(x-g)1 
^1 (x-0^ - /3^ r^ 




(24) 


Dorrance [3] used this form of the potential together with the slender body doublet 
strength to obtain stability derivatives for pointed bodies in supersonic flow. How- 
ever, as pointed out by Miles [26] , such a procedure is inconsistent and unneces- 
sarily complicated. Proceeding from equation (24) we will derive a more expe- 
dient form for the potential which will show the essential equivalence between 
Dorrance 's method and the Adams-Sears approach [20] and thus provide a con- 
sistent means of improving on slender body theory. 


For this purpose, we expand equation (24) for small distances from 
the body axis. Equation (24) can be separated into a steady part, 


(p^(x,r,e) 


cos 6 9 

2 7T 8 r 



0 


F(n 

^/ (x- -13^ r^ 




(25) 


which represents the potential of a body of revolution at small angle -of -attack, 
and an unsteady part. 


(x,r,e) 


cos 6 
2 7T 


8 
8 r 


x-^r 

/ 


0 


i U‘ jx-n ♦ F(g) 

7(x-^)^ - /3^ r^ 


d^ . 


(26) 


The angle -of -attack potential can be expanded by a slight modification of 
F. Keune's method for the zero angle -of -attack potential [27]. Integrating first 
the square root expression and taking the partial differentiation outside the inte- 
gral, we obtain for the zero angle-of-attack potential 


(p(x,r)=-f f F(^)-^ln [x- ^ + 'J(x-^f-/3^r^]d^ 

0 ^/(x-|)^ -/3^r^ 0 

8 I 

= - -r — r F( P) In [x-^ + N (x-4 ) ^ -/3^ r^ ] d| + F (x-;3r) * In /3r . 

d X 


8 



Expanding this expression for small r gives the familiar slender -body 
potential: 

X 

= F(x) In /3r - f F ( 4 ) In [2(x-4)] d| . 


( 27 ) 


0 


To obtain the next higher approximation, we repeat this process twice 
for the difference between the exact potential (pix^r) and its first approximation 

(x,r) ; thus, 


c^(x,r) - (x,r) = - / F(|)ln [x - 4 + -7 (x- ^ ^ 


x-/3r 


x 


+ F(x-| 8 r) • In jSr - F (x) In /3r + f F(^)ln[ 2 (x-|)]d^ 


x-/3r 


f F * ( 4 ) ~ < (x- 4) In [x-4 + (x-4) ^ r^ ] - J{x-^)^ - r^ > d 4 


r d 

+ F (x-jSr) In (3r - F(x) In /3r + J F*( 4 ) [ (x- 4 ) In 2(x-4) ~ (x- 4 ) ] d4 


x~l3v 


Bx 


J F’ (4 ) / (x- 4 ) In [x- 4 + n/(x-4)^ - r^] - n/(x- 4 )^ - r^ > d 4 


+ F’(x-jSr) /5r In /3r + F (x-/3r) In /3r - F (x) In /Sr 


X 


+ / F' ( ^) f (x-^) In 2 (x-^) - (x-^)] d| 


^ ^ jin fx-4 + 7(x-|)2 - r^] 


0X 


0 


J (x-4) ^/(x-4)^ - ) d 4 + F' (x-/3r) • /3r In j 8 r 


+ F (x-j 8 r) In jSr - F(x) In /5r 


X 


r F’( 4 ) i2LJJ 

8 x 0 x 2 


In 2(x-4) - J (x-4)^ d| 
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8 

8x 


- j (x-|) J (x-|)^ - /?* d I + F" (x-/3r ) . 


+ ^/ (x-i 


x-| + 'V(x-l)^ - /3^ 


3g^ 


In /3r 


+ F' (x-/3r) • |8r In jSr + F (x-/3r) In /3r - F (x) in /3r 


8x“ 


/ F’(4) 


ln2(x-4) - I (x-4)2 


d| 


(28) 


Using the following expansions , 


In [x-| + ^/(x-4)2 -/S^r^] = In [2(x-|)] - j -^^2 + 


and 

n/(x-^)^ - )3^ = (x-|) 


1- ^ ^ 
2 ( x -^)2 


we obtain from equation (28) the following second approximation 
(1) fi* r^ 8^ ^ 

0(x,r) - 0 Mx,r) = - — ^ / F’ (I) In [2 (x-|)] d | 


4 


F"(x) fln(jSr) - 1] . 


(29) 


Inserting this result into equation (25) gives, after differentiation with respect to 
r, the following near -field approximation for the angle -of -attack potential: 

i (x,r, 0) = F(x) + r [F" (x) ( In ^ - | ) 

R ^Trr 47 t ^ ^ 


dx 


j f F(|)ln(x-|)dn. 

0 


(30) 
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The unsteady part, equation (28) , likewise can be approximated by trans- 
forming it first into [28] 


^ cos e d 

0jU,r,6) = 


f 


x-/3r 


i u • (x-^) • F(^) 




cos 0 . „2 


2 x-/3r 


J F(4 ) cosh 


0 


'(s?) 


d| . 


This latter expression can be easily approximated for small r in the integrand 
and in the upper limit 


Therefore, the potential of the slowly oscillating body of revolution in 
supersonic flow, equation (24) , may be written in the following form: 

J(x,r,0)=j|^ + f(x,r;M)j cos 6, (31) 

where 

^ (x,r;M) = ^ [F"(x)^ln ^ - F"(0) In x 

- / F’"(|) In (x-^) d|] + r[F»(x) In ^ 

X 

- f F”(|) lu (x-4) d^ 1 . 

0 


x-j3r 


8x' 


f F ( I) cosh 


-1 


This result also could have been obtained by applying Fourier or Laplace 
transform techniques to the first-order potential equation and expanding the 
transformed solution for small distance from the body [7]. The present approach 
is more elementary and shows at the same time the essential equivalence between 
Dorrance's approach [3] and the Adams-Sears method [20]. 

Terms of orders like /3^ e^ln /3 e have been neglected in equation (31) so 
that it is consistent with the expressions given earlier for the pressure coefficient. 

For convenience we now define a reduced potential x (x,r) which follows 
from equations (4) and (31): 
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X (x,r) = 


6 • ip (x,r) + 6 • X (x,r) . 


(32) 


A 

cos 6 


The doublet distribution F(x) , as mentioned previously, is determined 
from the tangency condition at the body. Using equation (32) together with 
equation (7) , we may write the exact tangency condition as follows: 


dr 


-W(x) +R’(x) 


8x 


at r = R(x) , 


(33) 


where for a rigid body 


W(x) 


6 + 


x-a + 


Q'(x) ' 

2 7T 


(34) 


Q 

The product term in the tangency condition contributes O (e ) to the potential 
function and hence, for consistency, must be retained. 


When we apply equation (33) , an integral equation for the doublet distribu- 
tions F(x) is obtained. However, to try to solve this integral equation exactly is 
inconsistent with the order to which the potential is valid. The correct procedure 
is to apply the iterative method of Adams -Sears [20] . In keeping with the order 
to which the potential function, equation (31) , is valid, the Adams -Sears iteration 
consists of two terms only. Therefore, the reduced potential >; may be written 






(35) 


and likewise for the doublet distribution, 

F(x) = F^^^x) + F^^^ (x) . (36) 


In equations (35) and (36) , the first term is the slender-body value, while the 
second term represents a second -order correction due to Mach number and 
thickness effects. Substituting equations (35) and (36) into equations (31) and 
(32) and equating terms of like order gives 


(x,r) 

(x,r) 


2 7T r 


2 7T r 


^ f ^ \x,r, M) . 


(37) 


The body center-line motion function W (x) , following the Adams -Sears method, 
may be written 
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W(x) = 

(x) + (x) , 

\ 



where 



> 


(x) 

= 6 + (x-a) (5 


(38) 

(x) 

Q'(x) : 

- 2, ^ • 

) 



Notice that (x) represents the usual linearized rigid body motion, 

whereas 


(x) represents a higher-order term which enters through oscillation of the 
body -fixed coordinates and is of the same order as the product term in the tan- 
gency condition. Substitution of equations (38) and (35) into equation (33) and 
equating terms of like order gives the first- and second -order tangency relations: 


9r ' 

^ =-w' ^x) -hRMx)-^ 


) at r = R(x) 


(39) 


Explicit forms for the first - and second -order doublet distributions follow 
immediately upon application of the tangency relations, equation (39). 

(1) • " 
f' ’ (x) = 2Q(x) 6 + 2 (x-a) Q(x) 6 

O' 


F^^Nx) 


2 Q(x) Z^(x) 


>’^x) 

7T 


2Q(x) Z fx) - - (x-a) Q'^x)| 6 

1 7T 


where 


(40) 


^R<^^ = 4^ [Q"(x)(^ln +lj-2Q"(0) Inx 

X 

-2 J Q”'( 4 ) ln(x-|)d^ 

0 

((x-a) Q(x)]^ln ^ ^ + 1^ + 2 aQ"(0) inx 
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- 2 / --p- r(|-a)QU)] In (x-|) d| 




Q'(x) In + ^ - / Q’'(|) In (x-4) d| 


After some straightforward calculation using equations (40) and (37) , we arrive 
at the following expressions for the in-phase and out -of -phase cross flow pertur- 
bation potentials, valid to second -order: 


7T r 


1 + Zj^(x) 


in ^ ■ 


2 7 t" r 


2 


- r Yj^(x) 
A (x,r) 


(4i) 


^ [(x-a) +Zj(x)l - 


(x-a) Q' (x) 
2 7T^ r 


+ r In 


M . 


Vj(x) 


- r Yj(x) , 


(42) 


where 

, d^ m2 


Yr(x) = f Q"(x) + 2 Q"( 0) In X + 2 J Q"'( 4 ) In (x-^) d | ] 


■R' 


47T 


Yj (X) 


jL 

4i: \ dx 


2 f(x-a) Q(x)] -2aQ"(0) In x 


X 3 

+ 2 f —3 [(|-a)Q(|)] In (x-^) d^ 
0 ^ 

m 2 / 

+ j Q"(^^) In (x-^) d^ . 

7T J 
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STABILITY DERIVATIVE CALCULATIONS 


With reference to equations (11) and (12) , the four flutter derivatives are 


given by 


'dC 


N 


N 5 V 


6 ^0 


/ R C dx 

Q(l) 


(43) 


ac 


M 


yaa 


6 ^0 


■w) ( 


+ / R^R’C dx 

Q(l) J r» 


0 

'd c 


N 


N . 

6 


a 6 


6 — 0 


— y— f R c dx 
Q(l) ^ 


(44) 


(45) 


'a c 


M 


M 


a 6 




6 — 0 


/ R^R'C dx . 


Q(l) 


(46) 


The flutter derivatives are related to the stability derivatives by (cf, pp. 
16-17 of Reference 24) 


N, 


= C 


N 


a 


6 O' 


:: = c + c 

N. N. N 

6 a q 


M: M. M 

o O' q 


(47) 
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Numerous researchers in the past [3, 4, 5] have attempted to compute 
Mach number and thickness effects on stability derivatives by considering only 
the linearized form of the pressure coefficient: 


(l) + (2) 


C = - 2 !/) 

p 

1 


C = - 2 {X + lb) 

p X ^ 

2 


at r = R(x) . 


(l)+(2) 


(48) 


However, such a procedure is incomplete because the various quadratic terms 
in equations (18) and (19) together with the terms in the potential arising from 
the second-order tangency condition are all of comparable order to the thickness 
and Mach number terms in equation (48). 


To obtain all of the second -order terms we evaluate the pressure coef- 
ficient derivatives at r = R(x) as follows: 

0*'* - 2 (^ - -f)R’ W N 

1 ^ ^ 


c = - 2 ( A + + 2 [ f X^^^ + M^ il> + R(x) ] 

JP A A X 

2 


- 2 R'2 (x) 


M“ \ (1) M^ (1) 

1“) \ - T * -"'W 


(49) 


Needed in the evaluation of C and C are the slender-body potentials 

^2 


/'»andx"> 

• 

11 

IP 

2 7T 

ip^^\x,r) = 

Q(x) 

7T r 

II 

(x-a) Q(x) 
7T r 


2 27t 8x 


X 


/ Q*(x) In (x-l) d| 

0 


(50) 

(51) 

(52) 
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For conciseness we will evaluate the stability derivatives with the pitch- 
axis located at the nose of the body, i. e. , at a = 0. To compute the stability 
derivatives for any other pitch-axis location, we simply make use of the follow- 
ing axis transfer relations: 



where the subscript 0 denotes the derivatives calculated for the pitch-axis at 
the nose. The above transfer relations follow immediately from the definitions 
of , etc. , plus the following easily proved relation; 
a 

C (x,R) = C (x,R) -aC (x,R) . (54) 

p p P 

Before giving the stability derivatives, we first derive expressions for 
the rate of change with x of the in-phase and out -of -phase normal force coefficients. 
These are given by 


d C 




a 


d X 

d 


— ttR (X) Cp (x,R) 
1 


dx 

a 


^ j R(x) C (X,K). 

*0 


(55) 


(56) 


Appropriate substitution into equation (49) gives 
d C 


N 


a 


d X 


Q(l) 


d A^(x) 

d X 


2 B^(x) - 


Dr(x) 


(57) 
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d 

dx 


C + C 
N • N 
. a q. 


Aj(x) 


+ Aj^(x) -Bj(x) -Cj(x) 


-Dj(x) 


where , 

Aj^(x) = Q(x) < 2 + ^ Q" (X) In - 2 Q"(0) In x 


-2 / Q’”(|) In (x-l) d4 jj 
Bo(x) = ^ Q'(x) ^Q”(x) In ^ -Q”(0)ln x 

K 7T ^ 

X 1 


X 

-/ Q"'(4) In (x- 4 )d| 

0 J 


Dj^(x) = 2Q'(x) -Ij R'^x) + 


- Q"(x) 

7T 


and 

Aj (X) = Q(x)<^x+ txQ(x)]" In 


2 / U Q (4)]'" In (x-|) d| 

0 


.2 - QMx)ln™ 


+ n-2j Q”( 4 )ln (x-l) d| 


)> 


(x) = Q' (x)<[x Q (x)]" In 

1 ^ 7T 1 


5® Q (X) „ ^ 


-2 f [|Q( 4 )]”*ln (x-l) d| 
0 


+ — Q'(x) Q'(x) ( In +l)-2 f Q" (|)ln (x-|) d| 

7T V 4 7T / 


I I II III I I III II n nil II II wiiia 



Cj (X) = 

0 

- 2 Q" 

Dj (X) = 

0 


- fM2Q(x) + ^ xQ'(x)] [Q”(x) In 

7T ^ 

X 


(0) In X - 2 / Q’"(|) In (x-^) ] 

0 



i(xQ'’(x)] . 


Q (X) 
4 7T 


Q’(x) 


R^^x) 


In the evaluation of the second integral in the moment coefficient, i. e. , the 
moment of forces parallel to the body axis, the slender-body pressure coeffi- 
cient suffices. Thus, 

C (x,R) = - 4 R'(x) (59) 

C (x,R) = - 4 [xR(x)]’ . (60) 

^2 

0 


Performing the required integrations indicated in equations (43) through 
(46) with a = 0, we arrive at the following expressions for the stability deriva- 
tives for pointed -nosed bodies of arbitrary smooth profile: 

2 r 2 ^ 

=2 + Q’’(l) In - 2 f Q"’(x) In (1-x) dx 

JN 7T 4 7T 

a 0 


13^ 

7T Q ( 1 ) 



Q'(x) Q"(x) In 


/3^Q (x) 
4 7T 


dx 


1 i X 

- 2 Q"( 0 ) f Q’(x) In xdx - 2 / Jq'(x) Q'"(U In (x- 4 ) d^ dx 
0 0 0 




’(x) R’2 (x) dx 


7tQ(1) 


(61) 
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Q(l) 0 

0 




f Q(x) dx - ^ 

7T 


Q"(l) In 


g^(l) 


4 7T 


- 2 J Q*"(x) In (1-x) dx 




Q(l) 

1 X 


Q(x) Q"(x) In dx 

4 7T 


- 2 Q"(0) / Q(x) In X dx - 2 / / Q(x) Q’**(|) In (x-|) d| dx 

0 0 0 J 

+ — I 7 T. / X Q'(x) Q"(x) In dx 

irQ (1) [ •'q 4 7t 

1 lx 

- 2 Q"(0) / X Q'(x) In xdx-2jjx Q’(x) Q*"(4) In (x- 4 ) d^ dx 

0 0 0 


Q(l) 


1 - 




OfOI 


/ xQ'(x) RMx) dx •>■ 


7rQ(l) 


f Q'^(x) dx. 


(62) 




- 2 f [xQ(x)]”'ln (1-x) dx 


f Q(x) dx + / [2 Q' (1) + Q"(l)] In 

7T 1 4 7T 


4 M*' 


TT 


^ in * 1 - f Q"(x)ln(l-x) 

2 \ 4 7T / -' 


dx 


0 


11 


7tQ(1) 


1 2 1 
/ Q(x)Q"(x)ln dx - 2 Q”(0) / Q(x) In X dx 


lx 1 

- ^ f f Q'"(|) In (x-4) d| dx 

0 0 J 
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7tQ(1) 


1 X 


I i- f Q’(x) [xQ(x)]”ln dx 


- f /Q' (x) [4 Q(^) ] ’” In (x-^) dx 


0 0 
7tQ (1) 


1 x 




dx 


2 f f Q'(x) Q"(4) In (x-4) d^ dx 


0 0 


xQ7T)(/ [m’< 5(=‘) + ^ =‘'3'<x)j Q"(x)ln ^ ® dx 


-2Q"(0) / 


Q(x) + ^ X Q 


Mx)j 


In X dx 


1 x 

J 

0 0 




[M‘Q(x) + Y ^ Q 


'(x)j 


Q"'(^) In (x-^) d^ dx 


Q(l) 


/ 


IH- -) Q(x) - Mi - ^ X Q*(x) 


R’ (x) dx 


Q* (1) 

7T Q (1) 


(63) 


.c 

M • M 

) - -2- f ([2Q 

a 

'0 ' 

1 


-2 / fx 

Q (x) ] "' In ( 1-x) dx 

0 


4 [ 

o:m (^1, /3^Q(i] 

7T [ 

2 V 4 7T 


/ 3^Q (1) 

4 7T 


+ 1 ) - f Q" (x) In (1-x) dx 

4 7T y J 
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7tQ(1) 


/ f Q(x) [xQ X ]"ln dx 

\ 0 


1 X \ 

- 2 / /q(x) [|Q(|)]»"ln (x-I) d|dx) 
0 0 / 


M 


1 X 


f [(3(l)(ln //Q(x)Q"(«)ln (x- 5 )d{dx 




7tQ(1) 


f ^ 2 

f xQ(x)Q"(x)ln ^ (jjj 

I 0 


2 Q" (0) f X Q(x) In X dx 


1 x 


-2 f fx Q(x) Q"’(|) In (x-0 d| dx 


A 


0 0 
2 


IT Q (1) 


1 X 


1 2 
^ r X Q’(X) fx Q (x)] ” In P dx 

^ ^ 4 7T 


f /x Q'(X) [| Q ( 4 )] In (x-4) d| dx 


0 0 




7T Q (1) 


[/xQ'^x)(^ln + l^dx 


1 X 

-2 J JxQ'(x) Q”(|) In (x-^) d|dx 
0 0 


7 tq\i) I ‘1’^ M^Q(x)+ ^ xQ«(x)j xQ”(x)ln dx 


- 2Q"(0) / 


2 

M Q(x) + ^ X Q' (x) 


X In x dx 
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1 X 
2/ / 

o2 

Q(x) + 2 X Q’(X) 

0 0 


1 


7rQ(l) 

\ 2 ■ V / 

Q'^1) 

Q(l) 

7T Q(l) 

7T 


X Q”’ (4) In (x-4) dx 


Q'(x) + 2 7 t X R’ (x) 


dx 


(64) 


NUMERICAL RESULTS AND DISCUSSION 
Application to Specific Body Shapes 


To illustrate the effect of differing body geometry on the stability 
derivatives, numerical examples were worked out for three typical pointed 
bodies of revolution. These are as follows: 

1. Right-circular cone, R(x) = ex 

2. Convex parabolic ogive, R(x) = ex(2-x) 

€. X 

3. Concave parabolic ogive, R(x) = (1+x). 

These bodies were chosen, first, because they possess the three basic 
types of curvature, and second, to determine whether a similarity exists be- 
tween the geometric effects for bodies of revolution and low aspect ratio wings. In 
addition, other theoretical results are available for the cone [29] and convex 
ogive [7] with which to make comparisons. 

The actual evaluation of the expressions for stability derivatives given 
in the previous section is straightforward, although tedious because of the many 
integrals which occur. Even for the simple body profiles considered here the 
number of such integrals is very large, so that the main problem is one of book- 
keeping. Some simplification is to be had in the present case because all of the 
integrals reduce to four basic types, which are listed, together with their inte- 
gration formulas, in Table I. All required numerical values of the generalized 
formulas of Table I are listed in Table n. 

The final numerical formulas as functions of Mach number and body 
fineness ratio for the three bodies are given in the Appendix. The present 
theoiy^ can be applied to other smooth bodies, and analytical integration is pos- 
sible, provided the body area distribution Q (x) is a factorable polynomial. This 
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restriction arises in the treatment of integrals containing In Q(x). From a prac- 
tical standpoint, the labor involved would be considerable in treating any profile 
whose radius distribution R(x) is of higher degree than a quadratic. 


Discussion of Results 

Numerical results for the various stability derivatives according to the 
present theory were computed for the body shapes described above for three values 
of body fineness ratio, e = 0, 0. 05 and 0. 10. The case e = 0 corresponds to pure 
slender -body theory. A Mach number range from 1.2 to 3.0 was investigated, and 
a pitch axis location of a = 0. 4 was chosen. 


Plots of are presented in Figures 2 through 4, in Figures 5 and 

a a 

6, in Figures 7 through 9, and ^ in Figures 10 and 11. 

a q a q 

These plots show a significant influence of Mach number and body thickness on 
the four stability' derivatives. For the cone, excellent agreement with the "exact" 
first-order theorj'^ of Tobak-Wehrend [29] is obtained as shown in Figures 2 and 

7. In fact, it can be demonstrated by series expansion of equations (40) , (45) and 

(50) of Reference 29 that the present theory, equations (Al) - (A2) of the Appen- 
dix, results by retaining second -order terms only. No plots of or ^ 

a d 

C,, are given for the cone because the forces and moments are related to each 

q 

other by the simple exact relations [29] for a = 0 


C = - - (1 +e^) C 
M s' ^ N 

a a 

^ 0 


\ = -7(l + 0 


M . M 
a 


q/o 



(65) 

( 66 ) 


In Figures 3 through 6 the in-phase derivatives C and C for the two 

^0! ^ a 

ogives as given by quasi -slender body theory are compared with an exact method 
of characteristics solution [30] which superposes a small airgle -of -attack per- 
turbation on the nonlinear axisymmetric flow. The agreement for is good 

a 

up to a value of the hypersonic similarity parameter 2 M e of 0. 3. Fair agree- 
ment is obtained for the pitching moment coefficient slope . Better insight 

a 
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into the behavior of these derivatives is gained by comparing the local normal force 
distribution d Cj^ /dx of the two theories. These comparisons are given for vari- 
ous Mach numbers and fineness ratios in Figures i'J through 19. 


Referring to Figures 10 and 11 for pitch damping, we observe the following 
behavior. For a given set of governing parameters, the concave body, the cone 
and the convex body (in that order) are progressively more stable. Thus, in 
general, decreasing the body curvature increases the pitch damping, i.e. , 

C,, + C,, becomes more negative. This behavior is similar to the pitch 

a q 

damping characteristics as a function of leading edge cuiwature of a thin low 
aspect ratio wing in transonic flow as shown by Landahl [31] . 


Three nonlinear studies for the oscillating cone are given in References 
14, 15 and 32. Revell's analysis [14] is an extension to unsteady flow of the 
second -order slender-body theory of Lighthill. The analyses of Brong [15] 
and Hsu [32] , which are somewhat similar to one another, are based on a 
time -dependent perturbation about the nonlin^r Taylor -Maccoll solution [33] . 
The results of References 14, 15 and 32 for "*■ for a cone of 10- 

V “ 



degree semi-apex angle are compared in Figure 20 with those of quasi-slender 
body theory and the first-order "exact" results of Tobak-Wehrend [29] . The 
agreement between the present results and those of References 29, 15 and 32 
is remarkably good up to 2M e of about 0. 75, whereupon quasi-slender body 
theory begins to diverge sharply. The reversal of trends in quasi-slender body 
theory is similar to the behavior of Revell's second-order theory. By com- 
parison of the two with the "exact" first-order theory of Tobak-Wehrend, 
further support is given to Revell's conjecture that such trend reversals indi- 
cate the breakdown of second-order theory, and by inference, the breakdown 
of quasi-slender body theory also. 


The theory of Zartarian -Ashley for slowly oscillating bodies of revolu- 
tion in supersonic flow [7] , as mentioned earlier, also uses the Adams-Sears 
procedure to solve for the doublet potential. Their aerodynamic forces are 
obtained by applying 'linearized" momentum considerations, whereas in the 
present analysis these forces are obtained by integration of the pressures over 
the body surface. We observe that the present stability derivative results 
contain those of Reference 7 plus additional Mach number and body thickness 
terms. This disagreement, at least in part, is attributable to the use in 
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Reference 7 of the "linearized" momentum relation which, analogous to the use 
of the ’linearized" pressure coefficient, does not give all of the Mach number 
and body thickness terms. 


Finally, it should be noted that the present work is based upon an 
expansion of the cross-flow potential, equation (23), with respect to frequency. 
Therefore, the results imply, in addition to the limitations already imposed 
by linearized theory, that 

k M « (M - 1) 

which sets an upper limit to the frequency range and a lower limit to the Mach 
number range over which the results may be used. 


George C. Marshall Space Flight Center, 

National Aeronautics and Space Administration 

HimtsvUle, Alabama, November 24, 1965 
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FIGURE 2. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON 
NORMAL FORCE COEFFICIENT SLOPE FOR A CONE 
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FIGURE 3. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON NORMAL FORCE 
COEFFICIENT SLOPE FOR A CONVEX PARABOLIC OGIVE 




FIGURE 4. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON NORMAL FORCE 
COEFFICIENT SLOPE FOR A CONCAVE PARABOLIC OGIVE 
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FIGURE 5. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON PITCHING 
MOMENT COEFFICIENT SLOPE FOR A CONVEX PARABOLIC OGIVE 
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FIGURE 6. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON PITCHING 
MOMENT COEFFICIENT SLOPE FOR A CONCAVE PARABOLIC OGIVE 
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FIGURE 7. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON FIXED AXIS 
DAMPING IN PITCH NORMAL FORCE COEFFICIENT FOR A CONE 
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FIGURE 8. EFFECT OF MA-CH NUMBER AND THICKNESS RATIO ON FIXED AXIS 
DAMPING IN PITCH NORMAL FORCE COEFFICIENT FOR A CONVEX PARABOLIC OGIVE 
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FIGURE 9. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON FIXED AXIS DAMPING 
IN PITCH NORMAL FORCE COEFFICIENT FOR A CONCAVE PARABOLIC OGIVE 
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FIGURE 11. EFFECT OF MACH NUMBER AND THICKNESS RATIO ON FIXED AXIS DAMPING 
IN PITCH MOMENT COEFFICIENT FOR A CONCAVE PARABOLIC OGIVE 
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FIGURE 13. LOCAL NORMAL FORCE VARIATION FOR A CONVEX PARABOLIC OGIVE OF 
FINENESS RATIO 0. 05 AT A MACH NUMBER OF 3. 0 
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FIGURE 14. LOCAL NORMAL FORCE VARIATION FOR A CONVEX PARABOLIC OGIVE OF 
FINENESS RATIO 0. 10 AT A MACH NUMBER OF 1. 5 
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FIGURE 20. COMPARISON OF THEORETICAL RESULTS FOR FIXED AXIS DAMPING 
IN PITCH NORMAL FORCE COEFFICIENT FOR A 10' CONE 



Integral 

1 

/ m 

X In (ax + b) dx 

0 


1 

/ x“^ln(l -x) dx 
0 


f In X dx 
0 

1 X 

/ J x“ 4*\n (x-4) d ^ dx 
0 0 


I 


Integration Formula 
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TABLE II 
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APPENDIX 


Numerical Formulas for Stability 
Derivatives of Cones, Convex 
Parabolic Ogives and Concave 
Parabolic Ogives 


Cone ; R(x) = e x 

C_ = 2 - (2 + M^) 
N 

O ' 



Convex Parabolic Ogive; R(x) = ex (2-x) 


C 


N 


a 
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8. 0 In 
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J 
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01 
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/3e 


- 0.4190 M® + 3. 2762 


C,, +C^, \ = 3.0667 + 
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iSe 
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^ ' 13 € 
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cx 

Concave Parabolic Ogive: R(x) = — (1+x) 

A 


C„ = 2.0 - 8“’ (4.0 In - 6. 1874 ) + 1.5833 

r V J 
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Cjj ' - 1.4833 * 
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